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ABSTRACT 

Quasometry is precision measurement of celestial positions and apparent motion 
of very distant extragalactic objects, such as quasars, galactic nuclei, and QSOs. We 
use this term to identify a specific area of research, the methodology of which differs 
from that of general astrometry The main purpose of quasometry is to link the sub- 
milliarcsecond radio frame (ICRF) with the existing and emerging optical reference 
frames of similar accuracy, constructed by astrometric satellites. Some of the main 
difficulties in achieving this goal are discussed, e.g., the extended structures of quasar 
hosts, apparent motion on the sky optical variability galactic companions, faintness. 
Besides the strategic purpose, quasometry is undoubtedly useful for global astrometric 
surveys, as it helps to verify or even correct the resulting reference frames. There are 
two options of using measurements of distant quasars in a global astrometric solution: 
1) hard constraints embedded in the fabric of observational equations; 2) a posteriori 
fitting of zonal errors. The relative benefits and shortcoming of the two options are 
reviewed. A relatively small set of about 200 carefully selected reference quasars 
can go a long way in improving the astrometric value of a space mission, if they are 
sufficiently bright, stable, fairly uniformly distributed on the sky, and are defining 
sources in the ICRF. We present an ongoing program at the USNO to construct a 
quality set of optical quasars with the required properties and to enhance the ICRF 
with new sources in the areas where known, well-observed quasars are scarce. 

Subject headings: Astrometry - Reference Systems - Parallaxes - Proper Motions - Quasars: 
general - Methods: data analysis 
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1. Introduction 

The currently adopted International Ce- 
lestial Reference System (ICRS) is realized 
through the second release of the Interna- 
tional Celestial Reference FrameQ (ICRF2, 



Boboltz et aL]|20rO| ). The ICRF2 achieves a 
stability of the axes of ^ 10 yuas, and an in- 
dividual object accuracy down to ^ 40 yuas. 
To take advantage of this unprecedented ac- 
curacy, a future mission of optical astrome- 
try should directly measure the positions of 
ICRF2 sources along with reference and field 
stars. Depending on the expected intrinsic 
accuracy of the mission, these measurements 
can be used to improve the absolute accu- 
racy of the resulting optical reference system, 
or to verify the external accuracy with the 
independent data. In this paper, the first pos- 
sibility is discussed. The previous eff'ort in 
constructing an optical analog of the ICRF2 
preceding the anticipated global astromet- 
ric missions was directed toward compiling 
large catalogs of QSOs and AGNs ( [Souchay 
et al.||200"9l ), counting more than a hundred 
thousand objects, and the derivative opti- 
cal realization of the quasar-based reference 



frame (LQRF, [Andrei et aL||2009) , obtained 
by aligning the axes with the ICRS at the mil- 
liarcsecond level. Quasars are by no means 
ideal objects for optical astrometry. Various 
complications have to be dealt with, e.g., 1) 
optical variability, 2) extended structures, 3) 
apparent motion due to binarity, gravitational 
lensing, or jets. Even a small admixture of 
strongly perturbed sources can nullify all the 
perceived advantages of a quasar-based ref- 
erence frame, while the currently available 
information about the vast majority of LQRF 
objects is scarce. We describe an alternative 
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approach adopted at the USNO of selecting 
a much smaller set of a few hundred ICRF2 
quasars, carefully vetted and investigated at 
optical wavelengths. 

2. Global solutions with constraints 

The estimated performance of a global as- 
trometric mission is usually characterized by 
a mission-average error (e.g., of parallax), 
which is a very compact, but not so compre- 
hensive, characteristic of the propagating er- 
ror. A diff'erent approach is to consider the 
distribution of position-correlated error (of- 
ten called zonal error) on the sky, whether 
of accidental or systematic origin. Any dis- 
crete distribution of astrometric error on the 
sphere can be uniquely decomposed into a 
series of either vector spherical harmonics 
(for positions and proper motions) or scalar 
spherical harmonics (for parallax). Further- 
more, a global astrometric solution can be 
obtained in terms of the spherical harmonics 
rather than for individual objects. Due to the 
near-orthogonality of the spherical harmon- 
ics, such a solution is precise for any subset 
of the basis functions. Using this idea, the ac- 
curacy of a very large global solution can be 
quickly evaluated on regular computers, such 
as a standard laptop. This also allows us to 
perform a large number of mission simula- 
tions for diff'erent configurations of reference 
objects, single measurement precision, etc. 
Alternatively, since the eigenvectors of a dis- 
crete zonal error are the spherical harmonics, 
the covariances and the confidence of a given 
mission-average accuracy can be easily com- 
puted. The confidence levels are related to the 
uncertainty of the outcome, i.e., the sigma of 
the sigma. Following this recipe, the proba- 
bility distribution of the mission-average er- 
ror can be computed without multiple Monte 
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Fig. 1. — Distributions of mission-average paral- 
lax error of the SIM interferometer project with and 
without 23 ICRF quasars used as hard constraints in 
the global solution, but with otherwise the same in- 
put parameters. The remarkable improvement in the 
expected parallax accuracy and the tightening of the 
uncertainty of the mission outcome arises from the re- 
duced propagation of accidental zonal errors in low- 
order spherical harmonics. 

Carlo simulations of the entire mission. 

Figure [T] shows the distributions of the 
mission-average parallax error with 23 refer- 
ence quasars and without any quasars com- 
puted for the Space Interferometry Mission 
(SIM). Using only 23 quasars in the grid solu- 
tion as absolute zero-parallax constraints not 
only improves the mission expectancy, but 
practically eliminates the risk of a very bad 
outcome. Interestingly, doubling the num- 
ber of reference quasars results in a relatively 
modest improvement of the expectancy of ^ 
10%, and doubling again in an even smaller 
improvement. The reason for this weak de- 
pendence on the number of quasars is the 
"red" spectrum of accidental zonal error, i.e., 
that most of the error is carried by the low- 
order spherical harmonics. 

To fully realize this advantage, the quasar 
constraints should be hard-coded in the obser- 



Fig. 2. — Distribution of accidental parallax error 
(standard deviation) simulated for a global astromet- 
ric solution with about 80 grid quasars used as hard 
constraints in the condition equations. The quasars are 
shown with black dots. The error of unit weight is 
color-coded as shown in the color bar legend below 
the sky projection. 

vational equations. For example, the "hard" 
parallax constraints are implemented by re- 
moving the corresponding columns in the 
global design matrix. The benefits of such 
constraints are clearly seen in Fig. |2] where 
the standard deviation of the combined zonal 
error of parallax for a global astrometric mis- 
sion is depicted. The error expectancy is 
much smaller in the areas populated by refer- 
ence quasars than in the near-Galactic plane 
zone where ICRF quasars are not available. 
An alternative approach is to use the prior 
information (such as that the parallax of a 
quasar is practically zero) in a posteriori ad- 
justment of the solution. At the minimum, a 
constant zero-point of all parallaxes (which 
is the first spherical harmonic) can be deter- 
mined and subtracted. This approach is more 
cumbersome and less accurate. However, it 
may be safer if some of the reference quasars 
are astrometrically perturbed. 
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Fig. 3. — Archival HST image of the radio-loud 
quasar PKS 2349-014 at F606W. The map reveals a 
galaxy companion east of the optical core, separated 
by approximately 2 arcsec, and an extended asymmet- 
ric structure with tidal arms. 



PKS 2349-014 VLA map 8 cm 




Fig. 4. — VLA map of PKS 2349-014 reveals a per- 
fectly unresolved, single radio core. 



3. The dangers of quasometry 



Bahcall et al. (1995) investigated one of 



the nearest and brightest quasars, PKS 2349- 
014. Their images obtained with the Hubble 
Space telescope (HST) with the F606W filter 



(yellow-orange light) reveal a complex struc- 
ture consisting of several components (Fig. [3]) 
extending to at least 10 arcseconds. The unre- 
solved nucleus is the brightest component of 
this structure, while other components have 
estimated magnitudes between 17.9 and 21 
(Table [TJ. Even though the extended compo- 
nents are a few magnitudes fainter than the 
nucleus, they may strongly perturb the astro- 
metric determination because of their large 
extent and off'sets from the nucleus. Addi- 
tionally, the substrate galaxy is obviously off'- 
center with respect to the nucleus, and the 
faint companion galaxy is at a projected sepa- 
ration of ^ 2 arcsec. These eff'ects will cause 
a nontrivial photocenter shift. In particular, 
if we apply the center-of-mass (CoM) cen- 
troiding method to an object of this complex- 
ity, disastrous astrometric errors would en- 
sue, given in the last column of Table [TJ At 
the same time, the VLA map at 8 cm (Fig. 
|4]) confirms PKS 2349-014 as a high-quality, 
unresolved radio core without any extended 
features or companions. In this case, a bona 
fide reference quasar in the radio is a horrible 
source in the visible light. 

Using more distant quasars as reference 
objects for an optical frame will certainly al- 
leviate the problem, but only to some ex- 
tent. A quasar like the PKS 2349-014, but 
100 times more distant (which would be far- 
ther out than the most distant known quasar), 
would still be marginally perturbed for a mis- 
sion like Gaia. The situation is drastically 
better if a more sophisticated method of im- 
age centering than CoM is used. [Guyon et 
aT] ( |2006| ) performed a NIR adaptive optics 
imaging survey of nearby (z < 0.3) QSOs. 
Using the superb resolution of their images, 
they fitted the template point spread function 
(PSF) to all observed profiles. Estimating the 
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prospects of quasometry, the important result 
is that most of the images are quite close to 
the PSF in the core areas of approximately 1 
arcsec radius. As the majority of luminous 
QSOs are associated with either ellipticals or 
mergers, the host galaxies and other extended 
features become prominent, and eventually 
dominating, outside the cores. Therefore, us- 
ing specially adjusted PSF centroiding meth- 
ods restricted to the core signal should greatly 
mitigate the impact of extended features for 
many (but not all) reference quasars. 

[Camargo et al. | ( |2011| ) find a scatter 



of optical minus ICRF positions up to 80 
mas, which seems correlated with the X- 
band structure index. Their result implies 
that the radio loud quasars can not be used as 
reference objects in the optical wavelengths 
because of the perturbations in their photo- 
centers from bright extended features. The 
X-band VLBI images of most of the sources 
in their study show that the entire extent of 
the radio emission is typically less than 5-10 
mas. Assuming the optical and radio emis- 
sion are associated with the same active nu- 
cleus, from the astrophysical prospective, it 
is difficult to understand the nature of the 
large optical- VLBI off'sets reported in this pa- 
per. In our opinion, the most likely explana- 
tion for the large scatter is the local system- 



atic error of the UCAC3 catalog (Zacharias 



et al. 2010), which was not taken into ac- 



count in that paper. On a scale of 1 to 2 deg 
on the sky, systematic errors of proper mo- 
tions for stars in the 13.5 to 16 mag range are 
caused by the unmodelled systematic errors 
of the early photographic plates, which have 
scales of about 5° on a side. In particular, 
for the area north of -20°, UCAC3 rehes on 
Schmidt plate data, and systematic errors of 
proper motions up to 12 mas yr"^ have been 



found dRoeser et al.pOlO] ). Combined with 
the epoch diff'erence of about 8 years between 
the ESO, SOAR data from 2006-2008 and 
the UCAC mean epoch of 1998-2003, this 
gives rise to a systematic error of position of 
up to 100 mas, which will not decrease with 
the number of reference stars. For the area 
of the sky well covered by the SPM early 
epoch data ( [Girard et al.||2011 |), the system- 
atic errors in UCAC3 proper motions are es- 
timated to be much smaller, but the epoch 
diff'erence is larger (about 10 years), which 
could lead to 20 to 40 mas positional off'sets 
not taken into account in the paper by [Ca- 
margo et al. I ( |2011| ). Similarly large off'sets 



between the optical and radio positions of 



some ICRF quasars were found by ( [Zacharias 



& Zacharias 2008), where the possibility of 



large zonal errors in the optical reference sys- 
tem was discussed. 

An additional considerable complication 
arises from the optical variability of quasars. 
Taris et alT1 ( |2011| ) determined the light curves 



of 41 randomly selected QSOs in a single 
deep field. During more than 50 months of 
observations, most of the objects varied by 
up to 1 mag. The light curves in the G, R 
and / bands were strongly correlated, and for 
one of the investigated quasars, a correlation 
between the photometric variability and the 
astrometric location was detected. Such cor- 
related astrometric-photometric perturbations 
are to be expected if the photocenter of the 
image and the estimated centroid is aff'ected 
by the external structures, such as the host 
galaxy, because the extended components of 
the object are hardly variable. The situation 
is much better at radio wavelengths, where 
the extended features and inactive compan- 
ion galaxies are not visible. 'Fomalont et al. 



(2011) investigated four ICRF quasars and 
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Fig. 5. — Distribution of /-magnitudes of selected 
high-quality reference QSOs. 



found that some ICRF positions are domi- 
nated by a moving jet component and may be 
displaced up to 0.5 mas, reflecting the mo- 
tion of the jet. This probably sets the scale 
of the intrinsic apparent motion of the active 



nuclei. [Titov et aL] ( |2011| ) determined "proper 
motions", i.e., linear trajectories of 555 VLBI 
radio sources over some 20 years of observa- 
tion. Even though some of the quasars seem 
to have motions of up to 1 mas yr"\ the sys- 
tematic (correlated) part of the vector field is 
much smaller. Using their data, we redeter- 
mined the vector spherical harmonics up to 
order 3 for this set of sources, and set an up- 
per bound on the amplitude of the systematic 
vector field of 50 yuas None of the fitted 
spherical harmonics were statistically signifi- 
cant at the 3(T level. 

4. A precious set 

The best objects for quasometry are quasars 
that are both radio loud and optically bright. 
Furthermore, the source of radio emission, 
which mostly comes from the relativistic jet, 
should be compact and preferably without 
superluminal motion. The source of optical 



emission in the core, which, to a large extent, 
is generated by the accretion disk around the 
central black hole, should be compact and 
stable. The host galaxies should be faint 
and symmetric around the nucleus. These 
criteria seem to contradict each other. Ac- 
cording to [Sikora et al.| ( |2007| ), there exists 
a dichotomy of accretion luminosity versus 
radio loudness. They surmise that the black 
holes in ellipticals rotate faster, which leads 
to optically brighter nuclei. However, AGNs 
become radio quiet at very high rates of rota- 
tion. Presumably, high rates of rotation pre- 
clude coUimation of jets. This perhaps ac- 
counts for the sad tendency of most accurate 
ICRF sources to be optically faint. 

[Hamilton et al.| ( |2008D infer a "fundamen- 
tal plane" of QSOs on a sample of 70 objects 
with 0.06 < z < 0.46 observed with HST's 
WFPC2. This relation can be re- written as 



My(host) - My(nuc) = 



-0.697My(nuc) - 15.9 
(1) 



which includes the absolute V magnitudes of 
the nucleus and the host galaxy. According to 
this relation, the magnitude diff'erence is sta- 
tistically larger for more luminous nuclei, i.e., 
the intrinsically bright cores should be much 
less perturbed by the host galaxies. The host 
and the nucleus have equal magnitudes at 
My(nuc) = -22.8, but the diff'erence becomes 
2 mag at My(nuc) = -25.7. In this respect, 
the exceptionally luminous nearby quasar 3C 
273 at AMy = 3 is a good object for qua- 
sometry, but PKS 2349-014 at AMy = -0.2 
is bad. 

Many of the strategic purposes of optical 
quasometry can be fulfilled with a relatively 
small set of quasars of prime quality. For ex- 
ample, the alignment of the rotation and spin 
of an optical reference frame with the ICRS 
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Table 1 : Magnitudes and predicted astromet- 
ric offsets of optical components of PKS 
2349-014. The data in the second column 
are from ( [Bahcall et al.||1995| ), the data in the 
third column are derived by us. 



Component 




offset 




mag 


mas 


unresolved nucleus 


15.3 





wisps (tidal tails?) 


17.9 


150 


extended nebulosity 


18.0 


200 


companion galaxy 


21.0 


10 


host galaxy centered on nucleus 


18.6 






does not require a large number of common 
sources, but each source should be consistent 
at the optical and radio wavelengths to bet- 
ter than ~ 1 mas. We believe that selecting 
a limited sample of quality objects is a more 
promising approach than using as many un- 
vetted sources as possible. An ongoing ac- 
tivity at the USNO is addressed at vetting at 
least 200 QSOs that can serve as the foun- 
dation of the future radio-optical reference 
frame. A parallel activity is under way in Eu- 
rope dBourda et al.|2010D . 

We are using the following selection crite- 
ria: 

• Preferably high- accuracy ICFR2 source 

• Brighter than 7^17 mag 

• Photometrically stable, 61 < 0.5 

• Not a merger 
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• Not too close, z > 0.1 

• No nearby projections of stars or galax- 
ies 

• Nucleus is significantly brighter than 
the host galaxy. 

• Uniformly distributed on the sky 

The /-band magnitude distribution of the 
original selection of 518 ICRF and radio- 
mute QSOs is shown in Fig. |5j A second cut 
at this sample requires additional investiga- 
tion of individual objects involving archival 
data and new observations. In particular, the 
photometric campaign commenced in 2005 



from both hemispheres ( |Ojha et aL] |2009| ) 
goes on. It already resulted in the identifica- 
tion of highly variable objects, and of quasars 
with erroneous magnitudes in the catalogs. 

An extensive search for and perusal of 
images in the literature and in the digitized 
archives resulted in elimination of further 
several dozen objects. Some of the typical 
cases we have encountered in this process are 
depicted in Fig. |6} Some of the radio-loud 
quasars are projected too close to a bright 
star, others feature bright extended compo- 
nents or merger structures. Still, for many of 
the candidate reference QSOs, the available 
information is too scant for making a final 
decision. New good-seeing or AO images on 
large telescopes are needed for about a hun- 
dred of the original objects. 

Last but not least, the reference grid 
quasars should be as uniformly distributed 
on the sky as possible. The degree of "empti- 
ness" of the area around a given quasar is a 
significant figure of merit in our selection. To 
quantify this figure of merit, we developed 
a directional near neighbor statistic, which 




Fig. 7. — Distribution of defining ICRF quasars vet- 
ted for the future optical reference frame. Magnitudes: 
yellow = 12-14, blue = 14-15, green = 15-16, black = 
16-16.5, red = 16.5-17 mag. 



is more suitable for this task than the tradi- 
tional near-neighbor statistic (the distance to 
the nearest neighbor). If an object sits on 
the edge of a large void, such as seen in the 
lower left part of the sky projection in Fig. [7] 
around the Galactic plane, it should be given 
a high figure of merit even if another objects 
happens to lie close to it outside the void. 
The directional near neighbor (DNN) statis- 
tic is the mean of 12 angular distances to 
the nearest neighbors within 12 directional 
segments of equal width. More complicated 
rank orders can be constructed by combining 
the DNN statistics with magnitudes, since a 
strong preference is given to brighter quasars 
in our selection. 

We thank the USNO Editorial Board for 
helpful suggestions and a critical reading of 
the original version of the paper. 
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